On the onset of cosmological backreaction 
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Cosmological backreaction has been suggested as an explanation of dark energy and is heav- 
ily disputed since. We combine cosmological perturbation theory with Buchert's non-perturbative 
framework, calculate the relevant averaged observables up to second order in the comoving syn- 
chronous gauge, and discuss their gauge dependence. With the help of an integrability condition, 
the leading second order contributions follow from the first order calculation. 

We focus on the onset of cosmological backreaction, as a perturbative treatment is necessarily 
restricted to the era when the efi'ect is still small. We demonstrate that the leading contributions to 
all averaged physical observables are completely specified on the boundary of the averaged domain. 
For any finite domain these surface terms are nonzero in general and thus backreaction is for real. 
We map the backreaction effect on an effectively homogeneous and isotropic fiuid. The generic 
effective equation of state is not only time dependent, but also depends on scale. 

PACS numbers: 95.36.-l-x, 98.65.Dx, 98.80.-k, 98.80.Jk 



I. INTRODUCTION 

The accelerated expansion of the Universe has now 
been confirmed by various observations [H, [3, 0] ■ To un- 
derstand this mysterious phenomenon, different explana- 
tions have been suggested, e.g., dark energy in the form 
of a quintessence field or a modification of gravity on 
large scales. Among the many proposals, a cosmological 
constant (or vacuum energy) in the context of the infla- 
tionary A cold dark matter (ACDM) model seems to be 
the most attractive one, as it is simple and provides a 
good fit to all cosmological data currently available. 

However, the observational evidence is based on the as- 
sumption of a homogeneous and isotropic universe (the 
Friedmann-Lcmaitre- Robertson- Walker (FLRW) model). 
But, spatial homogeneity and isotropy are rather rough 
approximations for the Universe, only valid on scales 
larger than ~ 100 Mpc [1,01 • Thus, before assum- 
ing dark energy to be a component of the Universe, it is 
worthy to investigate the effects of local inhomogeneities 
and anisotropies. On the largest scales, the deviations 
from homogeneity and isotropy are tiny and a descrip- 
tion in terms of linear cosmological perturbations is very 
well justified. However, when the matter and metric fiuc- 
tuations enter the nonlinear regime (currently the largest 
observed non-linear objects extend to the 100 Mpc 
scale, among the most prominent the Shapely super- 
cluster and the Sloan great wall), the problem of aver- 
ages [l0( shows up, and as a consequence the observed 
cosmological parameters might actually differ from the 
actual parameters of the underlying cosmology. 

Let us pick the most important cosmological parame- 
ter, the Hubble expansion rate [6, 7] to discuss this issue 
in some detail. The idealized measurement of the Hub- 
ble expansion rate proceeds as follows. Take a set of N 



standard candles (in reality supernovae of type la) that 
sample a local (i.e., objects at redshift z <^ \) physi- 
cal volume V homogeneously. Measure their luminos- 
ity distances di (via magnitudes) and recession velocities 
Vi — czi and take the average [8| 



Ho 



1 ^ 
N ^ 



oo), this turns 



'Electronic address: li@physik.uni-bielefeld.de 
^Electronic address: dschwarz@physik.uni-bielefeld.de 



In the limit of a very big sample {N 
into a volume average 



In the second step, we neglect the effect of the light cone, 
but for z 1 the spatial average is a good approximation 
for an average over the past light cone, because the ex- 
pansion rate of the Universe is not changing significantly 
at time scales much shorter than the Hubble time. 

On the other hand, we have a theoretical object that we 
call the expansion rate, defined as = a/a, where a de- 
notes the scale factor of the background model. The issue 
in the averaging problem now is to establish the connec- 
tion between and . In linear theory, by construc- 
tion, the average Hq and the background H}^ agree if the 
volume V becomes large enough. However, due to the 
nonlinearity of the Einstein equations, cosmological per- 
turbations can affect the evolution of the average, which 
we often identify with the "background" Universe. This 
is the so-called backreaction mechanis ni |9|. [l(A [lU Il3 

Q E [H, [13, [H, [H, [M [M [13, [H, [H, [MlHia IH lli . 

As in the example of the Hubble constant, we typi- 
cally measure the value of a physical observable by taking 
an average of the observable in a domain of space-time. 
Other examples are number counts, correlation functions, 
power spectra, etc. Consequently, a comparison of the 
theory with observation needs to utilize averaged quan- 
tities. 

The influence of cosmological perturbations on the ex- 
pansion of the Universe shows up in many aspects, e.g., 
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for an underdense patch of the Universe the local expan- 
sion rate is naturally larger than its average. The task 
is to find the evolution equations for the averaged ob- 
servables and thus the effective equation of state of the 
Universe. 

If this effective equation of state WcS = Pcs/PcS < 
— 1/3, where pcff and p^g are the effective energy den- 
sity and pressure, the expansion of the averaged Uni- 
verse would accelerate. So, if backreaction (i.e., averag- 
ing) would give rise to a negative effective pressure, we 
might be able to explain the observed acceleration of the 
Hubble expansion. 

At the same time such a mechanism might be able to 
resolve the coincidence problem: why does the onset of 
acceleration happen around the present time? The an- 
swer of cosmological backreaction could be that originally 
tiny perturbations grow with time, and lead to the forma- 
tion of (weakly) nonlinear structures at the scales that we 
use to fix the cosmological parameters. Thus, any obser- 
vation that is based on physics in the (weakly or strongly) 
nonlinear regime might be influenced by the backreaction 
effect, especially the determination of Ho (based on local 
measurements), the Hubble diagrams from supernovae 
Type la (need data at low redshift and thus local infor- 
mation), the integrated Sachs- Wolfe effect, which can be 
mistaken for the nonlinear Rees-Sciama effect, and oth- 
ers. 

The discussion above does not imply that we have 
solved the puzzle of dark energy by the backreaction 
mechanism. It could well be that the backreaction effect 
is tiny and the nature of dark energy is indeed a constant 
vacuum energy density. Nevertheless, the backreaction 
effect is of interest, because with increasing experimen- 
tal precision, e.g., for the cosmic microwave background 
(CMB) effects of order 10~^ are measured [30|, and we 
must take the backreaction effect into account seriously. 

Recent research on the backreaction mechanism ex- 
plored two directions. One is to study the properties 
of the averaged physical quantities in the perturbed Uni- 
verse. In Buchert's work [H [H [13 (see [H, [H for 
a recent review), the averaged Einstein equations were 
derived in the synchronous coordinates with two fluctua- 
tion terms, the kinematical backreaction term (Q) u and 
the averaged spatial curvature {R)d- The behavior of 
the perturbed Universe thus depends on the properties of 
these averaged terms. In Buchert's approach, the average 
is taken over a physically comoving spatial volume. This 
seems to be the appropriate procedure for observations 
in our local neighborhood, such as for galaxy redshift 
surveys, which are currently limited to small redshifts 
{z < 1). However, for z > 1, spatial averaging does not 
reflect actual observations, as we observe the past light 
cone and not a spatial hypersurface. We therefore con- 
sider spatial averaging as a first reasonable approach to 
the cosmological backreaction problem. 

The second direction is to use cosmological perturba- 
tion theory to study the evolution of the perturbed Uni- 
verse, such as [l^,!!!!!!!!!!!!!!!!!!!!!!!^,!!! and 



references therein. All these works discussed the possibil- 
ity to explain the accelerated expansion of the Universe 
as the result of structure formation, without introducing 
dark energy into the Einstein equations. For example, 
in 20] , the Hubble expansion rate was calculated to sec- 
ond order in a dust (i.e., matter-dominated) Universe. 

Recent review and criticism on the backreaction mech- 
anism can be found in [2^, and references therein. 

Here, we synthesize these two lines of research. Be- 
cause doing perturbative calculations without averaging, 
we cannot obtain the global property of the Universe. 
Whereas, averaging without perturbative calculations, 
we cannot get quantitative information of the Universe. 
Therefore, in this paper, we calculate the averaged phys- 
ical quantities in the cosmological perturbation theory 
to second order, but without the need to use the metric 
perturbations of second order. In contrast to previous 
studies in the literature, we do not aim at calculating en- 
semble means or variances of spatial averages. It seems to 
us that the ensemble means of spatially averaged quanti- 
ties are of limited interest for the interpretation of actual 
observations. Our interest must be to quantify the back- 
reaction effect in the Milky Way's neighborhood (the lo- 
cal spatial domain of ^ (100 Mpc)'^, i.e., the domain 
used to measure Hq). We aim at predicting the amount 
of backreaction based on a measurement of the distribu- 
tion of matter density within that domain. As we show 
in this work, the knowledge of the peculiar gravitational 
potential on the boundary of a physically comoving do- 
main at some initial time allows us to predict the time 
evolution of the spatially averaged quantities (as long as 
the effect is small). 

Our paper is organized as follows. In Section 2, we 
introduce the concepts of the expansion, shear and rota- 
tion of the Universe, and use the ADM decomposition to 
rewrite the Einstein equations in terms of these quanti- 
ties. Following the averaging procedure of Buchert, we 
arrive at the averaged Einstein equations (Buchert equa- 
tions fT5|) for an irrotational dust universe and an inte- 
grability condition in Section 3. The integrability con- 
dition provides a consistency relation for the two back- 
reaction terms {Q)d and {R)d- With Section 4 we turn 
to cosmological perturbation theory, and solve for the 
first order metric perturbations \I' and x in the comov- 
ing synchronous gauge. In Sections 5 and 6, we calculate 
the backreaction terms {Q)d and (R)d, the averaged ex- 
pansion rate {9)d, the averaged energy density {p)d, the 
effective equation of state WcS and the square of the ef- 
fective speed of sound c^jj to first and second orders, re- 
spectively. Finally, we demonstrate in Section 7 that the 
backreaction terms {Q)d and {R)d, as well as and 
c^fj are gauge independent. Conclusions and discussions 
are provided in Section 8. 

In the following, the Greek indices run from to 3 and 
the Latin indices from 1 to 3, and we use units with c = 1. 
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II. KINEMATICS AND DYNAMICS OF THE 
EXPANDING UNIVERSE 

The standard FLRW model is based on the assump- 
tion of spatial homogeneity and isotropy of the Universe. 
However, these assumptions are not valid (even approxi- 
mately) at the scales on which structure formation hap- 
pens, i.e., the scales smaller with respect to the Hub- 
ble radius and sufficiently long after the matter-radiation 
equality. So, necessarily one must consider not only the 
expansion, but also the shear and rotation of the Universe 
in order to understand its kinematics thoroughly. 

A. Expansion, shear and rotation 

To describe the kinematics of the Universe, we need to 
calculate the gradient field of the 4- velocity = dx^/dr 
of comoving observers, where r is their proper time. We 
introduce the projection onto the spatial hypersurface de- 
fined by the comoving observers = (5^ -I- u^u^. Thus, 
the components of the gradient field of the 4- velocity de- 
fine the expansion tensor 

where 9 = ur\^, a^^ = h'^h'^iu^a-p] - \hapu^x)^ ^-^d 
Lu^y = h'^h^u\^a;p\ ^^c the expansion scalar, shear tensor 
and rotation tensor, respectively. 

In the following, we restrict our attention to an ir- 
rotational universe, i.e., = 0. Neglecting rotations 
seems to be a reasonable assumption in the context of 
infiationary cosmology, as there exist no seeds for vector 
perturbations and the conservation of angular momen- 
tum also implies that only nonlinear effects could lead to 
a generation of rotation. 

The metric of the inhomogeneous and anisotropic Uni- 
verse may be expressed in terms of synchronous coordi- 
nates 

ds^ = -dt^ + 5y (i,x)da;Ma;^ 

where t is the cosmic time, and x denotes the spatial co- 
ordinates. The corresponding nontrivial Christoffel sym- 
bols are 

r}fe = ^9'\93i,k + gikj -93k,i)- 

For an irrotational universe it is possible to use comov- 
ing synchronous coordinates (the observer is at rest with 
respect to the cosmic medium), i.e., = (1,0) and 
/i* — 6j. Thus, the nontrivial components of the shear 
tensor are 

<Jij = e,j - ^g,j9, (1) 



and we define the shear scalar as 

^'^l«-^(^}Ol-l<f)- (2) 

Furthermore, the nontrivial components of the expansion 
tensor become 

0^,^Tl^^g,,^o, = 0j = ^5%,,o. (3) 

From (l3|), we have 6 — J /J, where is the derivative 
with respect to the cosmic time, J = ^detg^ and detgij 
denotes the determinant of the metric. Thus, 

j = 0J. (4) 

B. ADM decomposition 

Having obtained all the quantities to describe the kine- 
matics of the expanding Universe, we turn to its dy- 
namics. For the dust Universe (in the comoving syn- 
chronous coordinates), the only nontrivial component of 
the energy-momentum tensor is Tq = — p, the energy 
density of dust. 

According to Arnowitt, Deser, and Misner i39|, the 
Einstein equations G^u = SttCT^^ in the present situa- 
tion can be decomposed into: 
the energy constraint 

R + -9]9l ^16ttGp, (5) 
the momentum constraint 

di-j = 0,^, (6) 
and the evolution equation 

0],o = -99)-R'j+At:Gp5), (7) 

where i?^- denotes the spatial Ricci tensor and R = R\ 
is the spatial Ricci scalar. Combining the trace of (|7l) 
with ((U and ([2]) leads to the Raychaudhuri equation [40] , 
which links the expansion and shear scalars together, 

9 = -\9'^ -2a'^ -A-KGp. (8) 

So far, we have not made any approximations apart 
from neglecting rotation and restricting the matter to 
dust. These equations are satisfied at any point in space- 
time. However, our observations do not allow us to mea- 
sure all the data that would be necessary to put a well 
posed Cauchy problem, but realistic observations deliver 
averaged quantities. In the next section we discuss the 
averaged properties of these equations, and in Sections 5 
and 6 we use cosmological perturbation theory to evalu- 
ate the averaged observables to first and second orders, 
respectively. 
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III. DYNAMICS OF FINITE DOMAINS 

In the last section, we set up the local dynamical equa- 
tions for a general irrotational dust universe, but realistic 
observations provide us with averaged quantities. We fol- 
low the averaging procedure by Buchert , and obtain 
the averaged equations of motion for the irrotational dust 
Universe. 

A particular choice in Buchert 's formalism is the focus 
on comoving domains. As long as we neglect the differ- 
ence between cold dark matter and baryons, the comov- 
ing synchronous coordinate system is uniquely defined. 
From the theoretical point of view, this appears to be the 
most natural set-up, given that we work in synchronous 
comoving coordinates. We argue in Sec. 7, that there 
are convincing reasons to pick this particular coordinate 
system, to be the one that is best adapted to the ac- 
tual observational situation. However, we should keep in 
mind that for observational purposes it is highly nontriv- 
ial to identify comoving domains. Nevertheless, we think 
that comoving domains are a reasonable approximation 
as the effect from peculiar motion is small compared to 
the Hubble expansion at scales beyond 100 h^^ Mpc. 



A. Averaging procedure 

The spatial average of an observable 0{t, x) in a phys- 
ically comoving domain D (with the dust particles) at a 
fixed time t is defined as [I4I 



(O)d^-^^ J(i,x)0(i,x)dx, 



(9) 



where Voit) = J{t, x)dx is the volume of the domain, 
and the boundary of the domain is assumed to be co- 
moving. Following this spatial averaging procedure, we 
calculate the averaged expansion rate {0)d as an exam- 
ple. 

From the definition of Voit), we may introduce an ef- 
fective scale factor an 



ap ^ f Vpjt) 
apo \ Vp„ 



1/3 



(10) 



where ap„ and Vp„ are the values oi ap and Vp at the 
present time. With the help of (|4]) and (jlOp we find the 
averaged expansion rate (6) p 



{0) 



p 



1 



(11) 



/ 0Jdx = 
Vp Jp Vp Jp Vp ap 

The effective Hubble expansion rate can thus be defined 



jdx=-^ = 3^. 
Vp 



Hr 



ap 
ap 



(12) 



The expansion rate {9) ^ and the Hubble expansion rate 
Hp naturally reduce to 3iJ and H = a/ a in the homo- 
geneous and isotropic case. 



An important consequence of the definition ([9]) is that 
the spatial average and the time derivative do not com- 
mute with each other. It is straightforward to prove a 
corresponding Lemma (commutation rule) ^l5] 



{O)o-{d)p = {oe)n-{o)p{e) 



p- 



(13) 



This Lemma is used to calculate the second order term 
of the averaged expansion rate {0)p in Section 6. 



B. Buchert equations 

With the definition of the spatial average ([9]) and the 
Lemma we yield the Buchert equations [T^j from 

averaging the Einstein equations (HI - and the Ray- 
chaudhuri equation ([5]), 



ap^ 




ap , 


/ a]j 


'dp 

ap 


A-kG , 
^ 3 



SttG 



-PcS, 



3poff) 



(14) 
(15) 



where p^g and p^s are the effective energy density and 
effective pressure of an isotropic fiuid, which read 



Poff = {P)P - 

1 

PcS 



1 



IGttG 



16nG 

{Q)p 



{Q)i 



{R)i 



_ 6kp \ 



{R) 



p 



6kp \ 



(16) 
(17) 



The expression {Q) p is the kinematical backreaction 
term. 



{Q)p EE -{{e-{9)pr)p 



(0)1) 



2{a^)p 
2{a^)p, 



(18) 



which consists of the variance of the averaged expansion 
rate and the averaged shear scalar. 

Equations (fT4|l and (fT5| express a highly nontrivial 
result! They closely resemble the Friedmann equations, 
but have been obtained without the assumption of homo- 
geneity and isotropy. What has been shown is that any 
irrotational dust Universe, averaged over comoving (spa- 
tial) domains appears to the observers to be a FLRW-like 
Universe. The "curvature term" kp /a'jj has been intro- 
duced to show that actually any FLRW geometry might 
be picked. Two geometries differ by their expressions for 
the effective fluid (flB]) and p7|) . Without loss of general- 
ity, we can thus take kjj — in the following calculations. 

This formulation of backreaction provides the link to 
arguments in favor of cosmological backreaction that 
have been put forward by one of the authors in [l^ . It 
has been argued that on the largest scales (where the 
cosmic principle applies) we can view the Universe as 
being described by a FLRW model filled with a single 
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isotropic, but imperfect fluid, i.e., we can then under- 
stand structure formation as an dissipative process that 
creates entropy and it has been shown in [l6| that the 
second law of thermodynamics imphes for an expanding 
dust Universe that Pcs < 0. 

From the Buchert equations, we see that the evolu- 
tion of the inhomogeneous and anisotropic Universe de- 
pends not only on the energy density, but also the back- 
reaction term (Q) d and the averaged spatial curvature 
{R) D- So it is quite important to know the values of 
{Q)d and {R)d- For instance, we find from (jlSp that 
if peff + 3peff < 0, i.e., {Q)d > '^'^G{p)d, the averaged 
expansion accelerates. In other words, the averaged Uni- 
verse can expand in an accelerating way in the dust era, 
even if the local expansion is decelerating everywhere in 
the Universe. Accelerated expansion of the averaged ex- 
pansion rate does not violate the strong energy condition! 

We calculate (Q) d and (i?) o in cosmological pertur- 
bation theory to both first and second orders in the next 
two sections. Furthermore, we can define the effective 
equation of state as 



From (|T4|) , (|T5l) and (122)) , we find the relation between 
{Q)d and {R)d (the integrability condition) [l^] 



(23) 



The integrability condition is an essential equation for 
the following calculations. Its advantage is that it can be 
used to any order in perturbative calculations, as it is an 
exact result. This is shown in Section 6, where we make 
use of the integrability condition to derive the second 
order terms of (i?) , {9)d and {p) d without using the 
metric perturbations of second order. If we had not made 
the choice fc_D = 0, we would need to replace {R) d in ((^ 
by {R)d — Bfeu/al), which would not change the solution 
of the integrability condition, as can be seen easily. 

So far, all our results are exact for an irrotational dust 
universe. In order to get quantitative information on the 
observed Universe, we turn to cosmological perturbation 
theory. We use the comoving synchronous gauge below. 



Wcff 



PcS 
PcS 



{R)d-3{Q)i 

mi 



IV. LINEARIZED EINSTEIN EQUATIONS IN 
(19) THE COMOVING SYNCHRONOUS GAUGE 



and the square of an effective speed of sound as 



Ccff 



PcS 
PcS 



(20) 



This effective speed of sound is the characteristic speed 
at which a small perturbation propagates through the 
effective fluid. An example would be a deformation of the 
boundary of the domain, or a perturbation caused by the 
introduction of some extra mass into the domain. Our 
effective speed of sound is different from the isentropic 
speed of sound. We calculate Wcs and c^g in Section 6. 



C. Integrability condition 

The Buchert equations contain two averaged quanti- 
ties, {Q)d and {R)d, which influence the evolution of 
the inhomogeneous and anisotropic Universe. However, 
these two terms are not independent, but can be linked 
by an integrability condition. 

In the irrotational dust universe, pressure is negligible, 
so the energy-momentum tensor is given by T^f = pu'^u^. 
From the covariant conservation of its time-like part we 
find the continuity equation 



In this section, we first introduce the metric perturba- 
tions and X and find the linearized Einstein equations. 
Solving these equations, we find the time dependence of 
^ and X- With the help of these solutions, we calcu- 
late {Q)d, {R)d, {d)D and {p)d to both first and second 
orders in the next two sections. 



A. Einstein equations for the perturbed Universe 

We start now from a spatially flat FLRW dust model 
with the scale factor a{t). In the synchronous gauge we 
write the first order linearly perturbed metric as [20| 



As' 



-di^ -I- a2(i)[(l - 2*)(5y + AjX]da;Mx^ (24) 



where 5* and x are the scalar metric perturbations at first 
order, Dij = didj — \5ij^ and A denotes the Laplace 
operator in a three-dimensional Euclidean space. The 
scale factor a in p4)) is not the same as the effective scale 
factor od defined in (jlOp . and their relation is shown in 
Section 6. 

From the line element (IM)) . we straightforwardly ob- 
tain the nontrivial Christoffcl symbols. 



P = 



(21) 



The space-like conservation is the Euler equation, which 
is trivial for an irrotational dust Universe, if expressed 
in the comoving synchronous coordinates. Taking the 
average of (PT|) and applying the Lemma (|13p we have 



D 



-{0)D{p)D^-'iHD{p)D- 



(22) 



aaSij — 2aa'^Sij 



pi — 



1 „ 1 

+ 2 



aaDijX 



D]dkX+^Dld,x^l:D,kd\, 



-DijX, 



(25) 



6 



and the components of the Einstein tensor 

2 



, a \ a ■ 2 / 1 . 



1 



-2d, ( vl' + -Ax 



2 



3* + 9-^- - (1' + ^Ax 
a \ 

1 .. 3a . 1 / 1 ■ 

2^+2^^^+^r+6^^ 



The energy-momentum tensor of dust becomes 

— P ~ P Pi 



(26) 



where p'"' and p^^^ are the energy density of the back- 
ground and at first order, respectively. 

We are now ready to obtain the hncarized equations of 
motion in the ADM decomposition. The different com- 
ponents at different orders are 
the energy constraint at zeroth order 



SttG 



,(0) 



and at first order 



3-* ^A ( * -f ^Ax ) - 47rGp(i)- 
a \ 6 



the momentum constraint 



^^\^ + -Ax 



0; 



the evolution equation at zeroth order 

2 



2^=0, 
a 



its diagonal (i = j) piece at first order 

3* + 9-^- - ^^A f * + I Ax) = 0, 
and its off-diagonal {i ^ j) piece at first order 



1 .. 3a . f / f ■ 



= 0. 



(27) 



(28) 



(29) 



(30) 



(31) 



(32) 



From the covariant energy-momentum conservation, 
we find at zeroth order 



and at first order 



+ 3^p(i) _ 3i,pio) ^ Q 
a 



(33) 



(34) 



Equation ((34l) has a first integral, 



C(x) 



3p(o) 



(35) 



which resembles the famous hypersurface-invariant vari- 
able (here for dust, expressed in the synchronous coordi- 
nates) 



C(i,x)^C-^Ax, 



(36) 



commonly used to characterize the primordial power 
spectrum ^41|] . 



B. Solutions for a, and x 

a. Solution for a. From (|33p. we have p^'^^'a^ = 

P^q'^o-q, where p^'^ and oq are the values of p*^"' and a 
at the present time. By means of ([27)). we find 



a 



2/3 



(37) 



We can see that a grows as t'^^^, which is the result of the 
spatially flat FLRW dust cosmology. But this does not 
mean that the perturbed Universe expands in the same 
way as the unperturbed one, because in the perturbed 
Universe, the meaningful scale factor is the effective scale 
factor an defined in (|10p . which, however, is not equiv- 
alent to a, and their relation is obtained in Section 6. 
So we cannot know the behavior of the expansion of the 
perturbed Universe in terms of the scale factor a. 

b. Solution for Wc first eliminate p^^^ from the 
equations of motion with the help of the first integral C 
from 



3«'*"a3 
^(,) ^ Sp^ ^ -^^^^ 



(38) 



This allows us to obtain an equation for 5". Namely, from 
|), dSIl) and (13HD, we have 



3* + 6-* = 47rGp' 



and using ()37|) . we obtain 



_ 127rGpfa| 



(VP + C(x)) , (39) 



•• 4 • 2 2 - 

3* + 7*-^* = ^C(x). 

So we find the solution for v]/ 

^'(x, t) = A(x)t2/3 + B(x)i-i - C(x), 



(40) 



where A(x) and i?(x) are constants of integration, i.e., 
functions of the spatial coordinates only. We can see from 
that vf consists of one growing mode Aipcjt"^/^, one 
decaying mode B{x)t~^, and one constant mode — C(x), 
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in which the free spatial functions must be fixed by the 
initial conditions. 

In the next sections, we see that only the time deriva- 
tives 5' and show up in the observables that are of 
interest to us. Also, we are only concerned with the evo- 
lutions of perturbations at late times, so we can neglect 
the decaying mode B{x.)t~^. Hence, only the growing 
mode yl(x)t^/'^ is of importance for the following calcu- 
lations. Thus, we get the time derivatives of 



We see from (|40p and ([45l) that both 4" and x grow 
as t^/"^ at late times. Because a cx t^^^, both ^ and 
X grow linearly as the scale factor a in the perturbed 
dust Universe. So if cosmological perturbation theory is 
valid, i.e., the perturbative terms ^' and x &re small, the 
scale factor should not be too large. In other words, a 
perturbative analysis (at any order) is restricted to the 
"linear" regime. 



-1/3 



-A(x)f 



-4/3 



(41) 



C. Relation between ^(x), C(x) and (^(x) 



c. 



Solution for x- From G* =0, we have 



3^ + 9_$ _ 
a 

3a 



-A 



2a' 



6 



Ax 



0. (42) 



Multiplying on taking the time derivative and 

inserting (P^ . we have 



3aa 



D]x 



2^f + 6-* 

a 



0(43) 



From ([371) and p4|l . the second part in ([43]) is 0, so we 

get 



D]X 



and thus 



D)x = qi^)t-'/' + E)i^)t 



where C'j(x) and £'j(x) are functions of spatial coordi- 
nates only. In the following calculation we neglect the 
decaying mode £'j(x)t~^, so 

D]x = q{^t-'/\ (44) 

From ([44]) . we obtain the solution for x 

x(x,t)=C(x)t2/3-fg(x), (45) 

where C(x) and g{K) are also the functions of spatial 
coordinates only. The function ^(x) does not carry phys- 
ical information, as it can be fixed by the residual spa- 
tial gauge transformation of the comoving synchronous 
gauge, and we utilize this freedom to set ^(x) = 0. Sub- 
stituting ([15]) into ([H]l . we can get the relation between 
C(x) and C](x) 



q(x) = -D]C{^). 



Finally, the time derivatives of x ^-re 



^C(x)i 



-1/3 



-4/3 



(46) 



We have obtained the solutions for 'J and x, but these 
two solutions are not independent due to ([29]) . From 
([^ . ([IT]) and ([H]), we have 

1 



i^A+-ACj =0. 

Because both A and C are functions of the spatial coor- 
dinates only, we have 



6 



(47) 



where K is a constant. In the spatially flat universe K = 
0, which is shown in the next section, and we find the 
relation between A and C, i.e., the relation between "if 
and X 

A = -IaC. (48) 
6 

Finally, with help of the evolution equations (|3ip and 
([5^ . we obtain 



''0 



(49) 



Let us note that at superhorizon scales C ~ C- As the 
amplitude of C at superhorizon scales is measured by cos- 
mic microwave background experiments, the magnitudes 
of the time independent functions C, C and A are thus 
fixed. 



D. Relation between C(x), (^(x) and the peculiar 
gravitational potential v^(x) 

The peculiar gravitational potential </5(x) is defined 
from the Poisson equation as [20| 



A^(x) = AttGp' 



(1)^2 



(50) 



From ([39]), using ([4T]) and ([48]), we have 

AC(x) = -127:Gp^^h'^/\ (51) 

So, with the help of ([51]) and ([37|), we obtain the 

relation between C(x), C(x) and (p{x.) 



'■0 



(52) 
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We see that the peculiar gravitational potential ((C(x) is 
just a linear function of the metric perturbation x- 

Therefore, we know the solutions for a, 4* and Xi which 
we use in the next two sections to calculate the first and 
second order contributions to the averaged physical quan- 
tities, and we focus on the investigation of the time de- 
pendence of the averaged observables. 



V. FIRST ORDER PERTURBATIONS 

In this section, we calculate {0)d, {R)d and {p)d hi 
the perturbed Universe, and these results are the first 
step to the derivation of the second order contributions. 
We do not calculate {Q)d, because it is a pure second 
order term, which is proven in the next section. 

For the first order calculations of averaged quantities, 
the integration measure J must be expanded to the first 
order as well, 



J = a^(l - 3*) =a^ (l- iAt^'^ 



1 



1 



as at late times, the decaying and constant modes are 
negligible. In the following, let us denote 



(O) 



Dl 



(53) 



which is defined to be the average on the background 
(i.e., J = a?). Watch out that the average is still over 
a physically comoving domain, which might have a dis- 
torted geometry, even on the background. Thus, for the 
first order calculations, the averages of the zeroth and 
first order quantities are 



Jdx 



Jdx 



(54) 



Therefore, the perturbation in J does not affect the first 
order calculations. 



A. Averaged expansion rate {9) d 

Substituting the perturbative connections in (|25p into 
([3]), we have 



0) = -5) ~ ^5) + \d\x, 



(55) 



and taking the trace, we get the perturbative expansion 
rate to first order. 



9 = 2,-- 3*, 
a 



(56) 



in which we have used the property = 0. Using ([37|) 
and (j4T|l , we obtain the averaged expansion rate {6) o as 
a function of cosmic time t and A, 



(0) -3" 3/vi/\ -2 



(57) 



From ([57|) . the first order perturbation decays as t~-^/^, 
which is slower than that of the zeroth order term 2/t. 
Therefore, the perturbation becomes more and more im- 
portant as the Universe evolves. However, this does not 
mean that this perturbation dominates at late times, as 
in the perturbative approach, we must restrict our anal- 
ysis to |(^')d| < 1. 



Averaged spatial curvature (R) d 



From (O and the trace of (O, we have 
R=l&nGp-e^ + 9)91, 



and thus 



R = UirGp - 9 



R = ~9^ - AO - 39)91 



(58) 



(59) 



By means of (155j), §7} and (gTj), we find to first 

order 



{R)d = 



40(A) 



Dl 



3i4/3 



(60) 



Closer inspection of ((60|) shows that 



1. {R)d has only a first order term while the zeroth 
order term vanishes, as the background metric is 
spatially flat. 

2. {R)d decays as t''^/^. From a oc t^/s^ go 
{R)d fx 1/a^ as the Universe expands [2^ . 

3. we may, with the help of (|T7)l . rewrite {R)d as 



{R)i 



40 



(—^{AC)di + Ky So we can see 



3f*/3 V 6 ' 

that the constant K in p7)) contributes to the av- 
eraged spatial curvature a term AOK/Sf^^^ oc 1/a^. 
We know that in the unperturbed k ^ universe, 
R = 6k /a^. Thus, 'iOK/3t'^/^ is just the back- 
ground spatial curvature term, and as we discuss 
the perturbations in the spatially flat Universe, 
this term vanish naturally. This is the reason for 
A = -AC/6 in 



C. Averaged energy density (p)o 

Similarly, from ([55]) . ([55]) and ([5^ . we have 



{())D+m) 

4ttG 



1 



(A) 



Dl 



6nGf^ 27rGi4/3- 



(61) 
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So we find for a domain overdense in average that (A) m 
is positive. At the same time, from ([57]) the averaged 
expansion rate is reduced. This is consistent with the 
intuitional understanding of the gravitational collapse, 
which decreases the expansion rate of the Universe. Also, 
from (j60|) we have a positive averaged spatial curvature 
for the overdense regions. 

We find in the first order perturbative calculations that 
only enters the expressions of {0)d, {R)d and {p)d, 
and the metric perturbation x does not show up. We 
show in the next section that cr^ = ^D^jX^iXi so x is re- 
lated to the shear of the perturbed Universe. This means 
that only the expansion influences the evolution of the 
averaged spatial curvature term and the averaged energy 
density in the perturbed Universe at first order. 

Let us finally note that the first order eontributions to 
{0)d, {R)d cind {p)d o,i"e surface terms actually, as 
we may write them as integrals of total derivatives 



{A)oi = 



6/^dx 



More surface terms show up below, when we turn to the 
second order perturbations. 



as [20|, |42, |43|, |4J| . However, in these papers, the metric 
perturbations of second order are always needed for cal- 
culations, and these calculations are always rather com- 
plicated and tedious. In this paper, we show how to 
obtain the leading terms of second order contributions to 
{Q)d, {R)d, {0)d and {p)d from the metric perturba- 
tions of first order only. 

We first prove that the kinematical backreaction term 
(Q) £) is a second order term, and then use the integra- 
bility condition, which is a crucial new input, to find the 
second order terms of {R)d, {0)d and {p)d- In these 
calculations, the shear and thus x show up in the ex- 
pressions. The effective equation of state Wcs and the 
square of the effective speed of sound c^g are also given 
to second order. 

Different from the first order contributions, at second 
order we have to consider the perturbation of the inte- 
gration measure J. Therefore, the averaged quantities of 
physical observables of different orders become 



('^^15 = (0«)m-3(0«vi/)^i+3(0«)m(*) 

= (O('^)m. (62) 



Dl, 



(2)\„ = 



D 



VI. SECOND ORDER PERTURBATIONS 

We move on to the second order perturbations of phys- 
ical quantities. Second order cosmological perturbation 
theory has been discussed widely in the literature, such 



We can see that at second order, the average of a first 
order quantity (0'"'^')d picks up two second order modifi- 
cations -3(0(i)^')di+3(0(^))di(^')di. In the following, 
we show that these modifications show up in the second 
order calculation of {6) £> . 



J 



A. Averaged kinematical backreaction term (Q) d 



Let us recall the kinematical backreaction term {Q)d defined in ([18] 



{Q)o^-mD-{e)l)-2{a^)o. 

Now we prove that {Q)d is a pure second order term. We show that the first part {9'^)d — and the second part 
{a^)D are both of second order. 

To calculate the variance {0'^)d — {(^)'d to second order, we write 



(63) 



where ^(o), 9'-^^ and 6l(2) are the zeroth, first and second order contributions to 9, respectively. 0'-°-' and 6*^^-' have been 
calculated in ((56|l . Using ([62|l . to second order we have 



^')d {9)1 - {{9^°^ + 0^'^ + 9^'^r)n {{9^"^ + 9^^^ + 9^''>)Df - {{9^'^f)D {9^'% - {{0^'^r)Di {9^''^}lim 



This means that the first part of (Q) d is of second order, however, for calculating it we do not need to know the 
detailed form of 9*^"^^ . All we need is 9 up to first order (see ((56)) ). 

Similarly, we calculate the average of the shear scalar {cr'^)D- From ((TJ) and ([55| we find at first order that 



"' "' 3 2 



(65) 



so (T* has no zeroth order contribution. Hence, using 



we have 



1 



18i2/3 



(66) 
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This means that {(t'^)d, the second part of (Q) d), is also of second order, and can be calculated by using the expression 
of X only at first order. 

So far, we have proved that both two parts of {Q)d are of second order, and consequently {Q)d is a second order 
term, but nevertheless can be calculated from the first order contributions to '5 and x- Using (pS)) . (|¥T1) and 

we get {Q)d at second order 

{Q)d = 6 ({^')m ml,) \{D)xDix)m = ^ [3 {{{ACr)m ~ {d^d.Cd^ d,C)ni) - 2{AC)l,] ^ ^,(67) 
where 

F ^ 1 [3 {{iAC)')Di - {d^d,Cd^d.C)Di) - 2{AC)l,] = ^ [3 {{d\d.C AC)) di - {d\d,Cd^d,C))Di) - 2{AC)l,] . 

F has only second order terms and is a function of spatial coordinates, and here we use (|48p to express the final result 
in terms of the variable C only. 

Some remarks on this result for (Q) d are in order. From (|67p we find that 

1- {Q)d, written in the form abbreviated -F, contains two second order terms, which are total derivatives and 
become surface terms when averaging. Meanwhile, the third term {AC)1^ is the square of a first order sur- 
face term, and thus its second order modifications in (|62p do not show up in F. Therefore, the kinematical 
backreaction is a function of the derivatives of C on the boundary of the averaged domain only. 

2. Because AC is a fluctuating term, and can be stochastically positive in some regions and negative in others, 
its average is expected to be negligible, if the averaged domains become larger (but are still on subhorizon 
scales). However, d^djCd^ diC and (AC)^ are positive definite, and therefore give nontrivial surface terms when 
averaging. Thus, {Q)d consists of these two surface terms on large scales. If they cancel, we can say that 
there is no kinematical backreaction at second order. In Newtonian limit, this cancelation was discussed in [i^ 
for periodic boundary conditions, and in [46l . [4^ for spherically symmetric spaces. In relativistic cosmological 
perturbation theory, this problem was treated in [isl. [48j. However, in general case, there is no reason for this 
cancelation. A review of this cancelation problem can be found in ,243 ■ 

3. We can see from (|67| that {Q)d decreases as t~^^^, which means that {Q)d oc 1/a. And we already know that 
(p^'^^)d oc and {R)d cx 1/a^. So the kinematical backreaction term {Q)d decays slower than {p)d and 
{R)d- Therefore, in the course of the expansion of the Universe the kinematical backreaction becomes more 
and more important in the effective energy density peff and effective pressure Peff- Of course, we should pay 
attention that {Q)d is a pure second order term, but {R)d has got a first order term, and {p)d even contains 
a zeroth order term, so we cannot conclude that {Q)d dominates the late time evolution of the Universe. The 
effect of the kinematical backreaction on the expansion of the perturbed Universe depends on its evolution with 
time t (or scale factor a) in its denominator, and also on the value of the surface terms in the numerator. 



B. Averaged spatial curvature {R)d 

In the last section, we have calculated the averaged spatial curvature {R)d to first order in ([SO]) , and in this 
subsection, we use the integrability condition ([23]) to get its second order part. But to do so, we first need to find 
the relation between the effective scale factor an, which is defined as the cubic root of the volume of integral in the 
averaging procedure in (|10p . and the scale factor a, which is defined to describe the expansion of the Universe in the 
perturbative metric (|24p , and then express ajj as the function of the perturbation A. Since the integrability condition 
is an exact relation to any order, and we have already got {Q)d to second order in (j67p . solving the integrability 
condition {a^^(Q)D) + o,%{aj-,{R)D) = 0, it is straightforward to obtain {R)d to second order. However, in the 
following, we show that we do not need to calculate ao to second order, but only first order is sufficient for our 
purpose. 

We may rewrite the integrability condition as 

6—{Q)d + {Q)d + 2—{R)d + {R)d = 0- (68) 
ao an 

Because {Q)d is already of second order, in the first term of we only need the zeroth order term of aj^jajj. In 
the third term, since (-R)d has no zeroth order term, we need the zeroth and first order terms of aD/oiD- 



11 



Prom (fTTj) and ((56)) . to first order we have 



so at first order 



aoita) 



ii) (70) 



Thus, if {A)oi is negative, the effective scale factor ajj grows faster than the ordinary result t^^^ in the unperturbed 
dust model. We can also see this from (pT|) that if (A) ui < 0, {p)d is reduced with respect to the background and the 
expansion rate gets positive modification, which is consistent with the intuition that the underdense regions expand 
faster than the overdense ones. 

Substituting (|69|) and (|67p into ((68|l . to second order we have 

^^^^ + 3 U " ) ^ ^Y¥7^-^- 

Solving this differential equation, wc find 

{R)d = D (^^y exp {2{A)n. {t'^' - t^)) + Ij^^, [l " -P {2{A)n, (t^^ - t^))] , (71) 

where D is the constant of integration, which is a function of spatial coordinates. For consistency, we must expand 
the solution up to second order 

^^)^ = ^ (^^0^' - ^iA)DiDtl + 5Ftl^') + ^ [2(A)DiDtt^' - 5F 

There is only one undetermined constant of integration D in the above expression. From (|60p . we know that {R)d 
has no zeroth order term, so D has only the first and second order terms, otherwise, the terms in the first bracket 
would give rise to a zeroth order contribution. Wc may write D = D^^' + D'-^\ where D'^^^ and D^^) ^rc the first and 
second order terms of D, respectively. Because {A)]ji is a first order term and _F is a second order one, {R)d can be 
written as 

- ^(^''^^0^' + D^X^' - mniD^'hl + hFtl") + -lj-^(2{A)niD^^hy' - 5F), (72) 

where the first term D^^Hq^"^ /t"^/^ represents the first order term of {R)d- It is matched with the first order expression 
PI| to fix 

i^(^) = (73) 

Substituting ([75)1 into ([7^ . and using (|^ . we find the averaged spatial curvature {R)d to second order, 

on ^^(2) c; 

= -^{AC)di + ^ + ^ [{{d\d,Cd^d.C))Di - {d\d,CAC))ni) + 2{AC)l,] , (74) 

where G*^) = D^^Hq^^ - 2{A) diD'^'^HI + 5Ftl^^. Neither D^^) nor G^^) can be fixed by matching to some known 
coefficients. However, the term G^^V^^^^ is unimportant at any time. Early on, —^{AC)di is a first order term, 
while G^^^ is a second order one, so it is negligible compared to the first term in ([7^ . Similarly, at late times, G^^^^ /t'^^'^ 
decays faster than -^^jj [{d^{djCd^ diC)) di — {d^ {diC AC)) di +2(AG)|);^], because both numerators are of second 
order, but the exponent of the denominator in G'^^ /f^^^ is the larger one. Thus, — gft/s (AC) di is the first order term 

of {R)d, which is the same as the resuh in and ^ [{d^{djCd^ diC)) di - (d' {d,C AC)) di + 2(AC)l,^] is the 
leading second order part at late times. Therefore, in the following calculations, we write {R)d as 

on K 

- -^{AC)di + ^ [{{d\d,Cd^d,C))Di - {d\d,CAC))Di) + 2{AC)l,] . (75) 

Thus, at second order, {R)d is again the function of surface terms. So with (|60|1. we find that {R)d is the function 
of surface terms at both first and second orders. 

In this subsection, we have extended the calculation of the averaged spatial curvature (R) d to second order by 
using the integrability condition. Its advantage is that we can do the second order calculation, without knowing the 
metric perturbations of second order. This is because the integrability condition is an exact result to any order, and 
we have got {Q)d to second order with only the first order perturbation theory. 
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C. Averaged expansion rate {9)d 

The second order perturbation of the expansion rate has been discussed in the hterature. For instance, in [20| . 
Kolb et al. used the metric perturbations of second order to calculate the averaged Hubble expansion rate and its 
variance. Here, the expansion rate 9 is defined in the same way as that in (20j . namely 9 = u\. However, the Hubble 

expansion rate in j2^ is defined as •\/87rG'(/9)i3/3, in contrast to the one in our work, which is defined in (fT2|) as 
Hd = = \{9)d. Below, the second order perturbation of the expansion rate is calculated, but without using the 
metric perturbations of second order again. We also show that our calculation is consistent with the result in [20| . 
From (|55|) . we have 

{R)d = -{9^)d - mo - i{9]9l)D = -2{9^)d - mo - 6{<j^)d- (76) 

Since {R)d has already been calculated to second order in ([75]) . {<t^)d is a pure second order term, and we know the 
zeroth and first order terms of {9)d from ((57)) . we can obtain the second order perturbation of (6')£) from ((76|) . 
Using we expand 9 as 

9 = 0(0) + + 0(2) ^ 3^ _ gvi/ + (?(2) ^ 2 _ 2A ^(3)^ 

a t ti/3 ^ ' 

so to second order 

02-4 8A AA^ 40(2) . 2 2 A -^^^ 

^ -^-^ + ^^ + ^' ^--^+3^+^^^- (78) 
Substituting ^ into and using dHBl) and (gHl), we find 

^=^-40(2)-f0(2)_M!-J,,;.^Z.^-^. (79) 

We see from ([7^ that R has both the first and second order terms, so at second order, the first order term gives 
two additional second order modifications when averaging as shown in (|62p . Therefore, the average of the spatial 
curvature R to second order is 

on o 1 

= -^(AC)z,i - 4(0(2))^, - - ^ [ll((AC)2)^i - 10((AC))2,, + ?,{d^d,Cd=d,C)Di] .(80) 

Above, the Lemma allows us to write {9^'^^)di = {0^'^^)'di second order. Matching with ([73]) yields 

(^(2))^^ + + ^ [3((AC)2)^i +4(a^9,C9-'a.C)z,i] = 0. 

Solving this differential equation provides us with the second order contribution to the averaged expansion rate (0)_d 

{e^^^)Di = [m^Cf)Di+A{d^d,Cd^d.C)Di] + |, (81) 

where / is the constant of integration, and at late times the term I /t^ is negligible without doubt. Therefore, we find 
the averaged expansion rate (0) d to second order, 

{9)d = (e(")>Z5 + (0(l))D + (0('))D 

= ^ + ^(AC)z.i - ^ [4 {{d\d,Cd^d.,C))oi - {d\d,CAC))ni) + 7{{AC))l,] , (82) 

where the first order term 0(^) also contributes via averaging to the second order result (see (|62I) V Straightforwardly, 
from ((22|) . the averaged Hubble expansion rate is 



13 



We can see that (9) d and (H) o are also functions of surface terms at both first and second orders. 

Finally, we show that the result in ((82| can also be obtained by using the metric perturbations of second order 
in m, 



ds^ = a2(^) 



drf 



1 - 2*(i) - A 



+ A, (x« + \x^'^) + \ (9.xf ^ + 5,xr + x!f )) d.M.^- 



, (83) 



where r\ is the conformal time, Vl/'-^^ and x'"'^'' ai'e the first order scalar perturbations (the same as in this paper), ^'f^' 



,(2) 



and x^^^ are the second order scalar perturbations, xl is the second order transverse vector perturbation, and 

is the second order transverse and traceless tensor perturbation. Ignoring decaying modes, for a dust Universe, 
y^-^ and ^'(■^^ are also given in [20j . 



*W(77,x) = ^^(x) + ^A<^(x), x(')(r/,x) = -^ 
5 18 c 

*(2)(^,X) = -|^2(^)_|^5.^(^)5^^(^)^^ 



-<y5(x), 



10, 



(A^(x))^ - -a'a,^(x)a^a.(^(x) 



(84) 



From 



and 



in the comoving synchronous gauge, we have 
1, 



a a 



a 



2a 



2a 



where ' denotes the derivative with respect to the conformal time 77, and in the dust Universe, t = {ai^-qY / {27t^) . So 
using (|52p. we can straightforwardly get the average of the expansion rate 0, 



bal 



t 3ii/3 
^1/3 



54in/^ii/3 



[(9*(Ca,C))z3i + (CAC)z5i - 6(C)z5i(AC) 



42 



4 ((9*(5,C9^5,C))i3i - {d\d,CAC))Di) + 7((AC))?,i] 



(85) 



Thus, we find that the leading second order term in (|85p . which we get by using the metric perturbations of second 
order in [2^, is the same as that in ([82)) . One can see as already argued for the case of {R)d that the subleading 
second order contributions show the same time dependence as the first order term. Thus, it is justified to neglect the 
subleading terms as they can never (in the perturbative regime) overcome the first order ones. 



D. Averaged energy density {p)d 



Similarly, from pT|) and (I14p . we have 



{P)i 



{Q)d + {R)i 
IGttG 



Using (p7)l . ([75)1 and ([5^ . we get the averaged energy density to second order, 
1 



{p)i 



1 - ^(AC)i3i + ^ (2 {{d'{djCd^d,C))Di - {{d\d,CAC))Di) + 7{AC)l,) 

Z Zo 



(86) 



(87) 



and {p)d is a function of surface terms at both first and second orders too. 

Up to this point, we have obtained all the averaged quantities {Q)d, {R)d, and {p)d to second order, and we 
only have to use the first order perturbation terms \E' and x (without the necessity of knowing the metric perturbations 
of second order). These simplifications are based on the integrability condition and the fact that (Q) u is a pure second 
order term. Unfortunately, we cannot extend this method to higher orders. For example, if we go to third order, to 
calculate cr^, we would need 0^ = Fg^ to second order, and thus we must know the metric perturbations of second 
order. An exception is the irrotational and shearless universe, which is merely the FLRW model. 
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E. Effective equation of state and square of effective speed of sound 



a. Effective equation of state Weff ■ From (fin)) , using ((^7|) . ([7^ and ([57)) . we obtain the effective equation of state 
to second order, 



18 



t^/^ (88) 



Equation ((88)) is the perturbative expansion of the effective equation of state, which is both time and domain depen- 
dent. We can find its time dependence at different orders, and domain dependence in different coefficients. From (j69p . 
it is easy to get t as the function of the effective scale factor od, and at late times (i.e., < to) w^n can be rewritten 
as 



Wcff 



'(2)^2 



5{AC)di tl'^ 
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-an 



1 +4/3 



"0 
On 



(89) 



Equation (j89p is also a perturbative expansion of the effective equation of state, but in terms of au, which is of more 
interest than a. We can see that the second order coefficients in (|88p and (|89| arc different, this is because an is not 
proportional to t'^^^, and thus the second order coefficient picks up nontrivial contributions from the first order one. 

Therefore, WeS vanishes at zeroth order. This is different from the cosmological constant, for which wa = —1. 
Consequently, in a perturbative framework, the backreaction mechanism cannot induce accelerated expansion of the 
Universe as that would imply WcS < —1/3. Nevertheless, the cosmological perturbations allow us to investigate a 
possible change of the expansion rate of the averaged Universe that might in the later nonlinear stage lead to the 
accelerated expansion of the Universe. We discuss this on both small and large scales. 

Using ([5T|) . we may rewrite WeS as 
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Firstly, on small scales, {p^^^)di may significantly deviate from 0, so the first order term dominates the value of 



WeS, l-e., WcS 



(p(i))£)]^t2 gg^jj ggg fj-om (|90| that if {p''^^)di < 0, which means that the energy density is 
underdense locally, Wcs is negative, and since WcS oc t^, this effect will be more and more influential as time goes 
on, and might cause the accelerated (averaged) expansion of the inhomogeneous and anisotropic Universe. Of course, 
with the above expression we can trace the evolution only for small perturbations. Once they arc in the nonlinear 
regime, our approach fails. 

Secondly, for large averaged domains (but are still on subhorizon scales), like the discussion on F, the average of 
AC is expected to become negligible, since it is a fluctuating term, and only the two surface terms give nontrivial 
contributions. Therefore the value of WcS is dominated by these two second order terms on large scales. 
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We can see from that the sign of WcS depends on the contrast of the two surface terms. It vanishes for certain 
boundary conditions, see |23.145|. l46l. [i?!. [Isj . However, we think that these boundary conditions are not natural 
and that the generic case for a finite domain in the Universe is that the effective equation of state is given by a finite 
surface term, that might be positive or negative, depending on the details of the fiuctuations on the boundaries. 

An important lesson that we learn here is that the cosmological backreaction introduces an effective equation of 
state, which is not only time dependent, but also scale dependent. 



b. Square of the effective speed of sound c^^. Similarly, for the square of the effective speed of sound, we have 
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Like WcS, on small scales. 
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So if the cosmic medium is overdense locally, c^^ > 0. But we also see that c^g can be negative in underdense regions. 
Usually this suggest that some damping is going on, which is related to dissipative phenomena and the increase of 
entropy. These aspects will be investigated in more detail elsewhere. 
On large scales, the second order terms dominate and we find 

els = ^ {{d\d,Cd^d^,C))Di {d\d.,CAC))Di) t^l^. 

Also the sign of the square of the effective speed of sound depends on the contrast of the two surface terms. 

To summarize this section, we can see that all studied physical quantities, {Q)d, {R)d, {(^)d, Ho, PcS, WcS and 
c^g, can be expressed as functions of surface terms at both first and second orders. Thus, to know the values of these 
physical quantities, we do not need to know anything about the interior of the averaged domain. Only the physical 
information (i.e., the values of C, namely the peculiar gravitational potential ip, and its derivatives) encoded on the 
boundary of the domain matters. 



VII. GAUGE DEPENDENCE OF THE 
AVERAGED QUANTITIES 

We should finally discuss the gauge dependence of the 
averaged physical observables. In [43,1431, the gauge in- 
variance of physical observables at different orders are 
discussed in detail: to second order, a quantity is gauge 
dependent unless its zeroth and first order terms van- 
ish, apart from the trivial cases that it is a constant 
scalar field, or a linear combination of products of Kro- 
necker deltas with constant coefficients on the back- 
ground. Thus, we know that (Q) is a gauge invari- 
ant quantity, since it has only the second order term. 
Secondly, {R)d, WeS and c^g, which have the first order 
terms, and {0)d, Hjj and pos, which have both the ze- 
roth and first order terms depend on the gauge choice. 
However, the first order terms of {R)d, w^eff and c^fj are 
gauge independent as well. So in summary, we conclude 
that all leading terms of all physical observables are gauge 
invariant, while the higher order ones are not. 

This raises the question of the coordinate system. It 
seems to us that the comoving synchronous gauge is 
very close to the coordinate system of a real observer. 
Real astronomers and their telescopes are comoving with 
matter (we neglect the difference between baryonic and 
dark matter here), they use their own proper time in 
all of their observations and regard space to be time- 
orthogonal, which defines precisely the slicing and gauge 
that we use throughout this work. 



VIII. CONCLUSIONS AND DISCUSSIONS 

In this paper, we use both the Buchert equations and 
cosmological perturbation theory to study the evolution 
of the perturbed dust Universe in the comoving syn- 
chronous gauge. We investigate the possibility to explain 
the accelerated expansion of the Universe without dark 
energy. We calculate the averaged kinematical backreac- 
tion term {Q)d and the averaged spatial curvature {R)d, 
and find that (Q) jj is a pure second order term, and 
{R)d has both the first and second order terms. As we 
use a perturbative approach, these terms can only affect 



the evolution of the Universe perturbatively, and thus 
we can only hope to find an onset of the cosmological 
backreaction mechanism in this work. In some circum- 
stances, for example, on small scales, on which the effects 
of fluctuations of the energy density (p) d are significantly 
nonzero, the backreaction mechanism should not be ne- 
glected carelessly. 

We conclude that cosmological backreaction is for real 
and that it can both increase or decrease the expansion 
of the averaged Universe, depending on the averaged do- 
main under consideration. Thus we argue that the effec- 
tive equation of state of the Universe is time and scale 
dependent and so is the square of the effective speed of 
sound. 

We find in our perturbative approach, that all physical 
quantities are surface terms or squares of surface terms. 
This suggests the conjecture that also a nonlinear treat- 
ment would find only functions of surface terms. 

Another point of our paper is that we show in Section 6 
how to calculate the averaged quantities to second order 
for the leading growing mode, but use only the metric 
perturbations of first order. This is a consequence of 
the integrability condition, which is valid to any order. 
And this greatly simplifies the perturbative calculations, 
which are usually done with help of the metric perturba- 
tions of second order in the previous papers. 

Finally, we discuss some observational aspects of the 
perturbative calculations in our paper. The problem is 
how to measure the perturbations, such as 5* or A. 

The fluctuation amplitude is defined as 




Here we only calculate cr|, to second order, so the terms 
containing p*^^) are negligible. (Usually this amplitude is 
defined at the distance of 8 h^^ Mpc, but we can also 
define it for any domain D.) From (|6ip . we have 
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So by measuring cr|,, we can know the value of 

Secondly, a measurement of the averaged energy den- 
sity {p)d in ([57]) . the averaged expansion rate {0)d in 
([82)1 and the averaged spatial curvature {R)d in ((75)) . 
together with the measurement of a^, would allow us 
to find {d'^djCd^ diC) Di- Because all of them are physi- 
cal observables (at least in principle), from these mea- 
surements we can know the values of ((AC)^}^)! and 
{d^djCd^diC)£,i. And since the other averaged physi- 
cal quantities {Q)d, w^s and c^g are the functions of the 
former, we can obtain the quantitative information of all 
these averaged terms. This will help us to understand the 
evolution of the inhomogeneous and anisotropic Universe 
and its relation to what we call the background model on 



a much deeper level. It also demonstrates that we do not 
need to retreat to a statistical treatment of the back- 
reaction effect, but we can try to design an experiment 
to directly measure its sign and magnitude in the Milky 
Way's neighborhood (e.g., the local ~ (100 Mpc)^ 
domain) . 
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